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The effect of string and quantum gravity inspired minimum-length deformed quantization on a 
free, massless scalar field is studied on de Sitter background at the level of second quantization. 
Analytic solution of a field operator is obtained to the first order in deformation parameter. Using 
this solution we then estimate the two point and four point correlation functions (with respect 
to the Bunch-Davies vacuum). The field operator shows up a non- linear dependence on creation 
and annihilation operators, therefore the perturbation spectrum proves to be non-Gaussian. The 
correction to the power spectrum is of the same order as obtained previously in a similar study that 
incorporates the minimum-length deformed momentum operator into the first quantization picture 
^SJ . and then proceeds in the standard way for second quantization. The non-Gaussianity comes at 

the level of four point correlation function; its magnitude appears to be suppressed by the factor 
' ^ exp(— 6A'^), where A'^ is the number of e-foldings. 
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I. INTRODUCTION 

A striking consequence of inflationary picture is that all structures we observe today owe their existence to the 
very small quantum fluctuations (of Planck length and even smaller) occurring during the inflationary epoch. So, 
it represents a very instructive example of the strong connection between Planck scale physics and the large-scale 



I i universe [ij, I4I . In view of the high precision cosmological measurements the study of possible imprints of quantum 
gravity on the cosmological perturbations could open an experimental window for testing of Planck scale physics. The 
aim of this paper is to study further the influence of minimum-length deformed quantum mechanics 

^ : [Q, p] = *(i + f^p') , (1) 



upon the inflaton perturbation spectrum. (We will assume natural system of units h ^ c ^ 1 throughout this 
paper.) In previous studies pertaining to the cosmological perturbations most of the results were obtained 



\l ' by incorporating the minimum-length deformed momentum operator in the first quantization picture and then the 
standard scheme for a second quantization was applied. To date, in the framework of this approach a substantial 
amount of theoretical work has been devoted to various phenomenological questions. The following obvious step is to 
complete this study by undertaking the minimum-length deformed quantization at the level of second quantization. 
^ , A Specific objective of this paper is the quantization a massless scalar field on the de Sitter background with respect 
Jv>( to the minimum- length deformed formalism, Eq.([T]). For this purpose we shall apply the theoretical framework for 
[ second quantization with respect to the minimum- length deformed quantum mechanics put forward in fldi [ll| . First 
a • we write down the corrected Hamiltonian, then solve Heisenberg equation of motion for field operator to the first order 
in /3, and then estimate two point and four point correlation functions with respect to the Bunch-Davies vacuum. 

Let us now somewhat quantify our discussion. The minimum- length deformed quantum mechanics Eq.([T]) underlies 
the generalized uncertainty relation 



5QSP > ^ (1 + PSP') , (2) 

proposed originally in the context of perturbative string theory as a consequence of the fact that strings cannot probe 
distances below the string scale (string length) [12]. This relation immediately imposes lower bound on position 
uncertainty SQmin = y/P- Further, Eq.(l2]) was discussed by combining the basic principles of quantum theory and 
general relativity in the framework of Gedankenexperiments 13]. The parameter /3 is set by the quantum gravity 
scale (either by string length or Planck length which are of the same order) ^/]3 ^ Zp ~ 10^"^'^ cm. The deformed 
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position and momentum operators in Eq.(IT]) can be solved in terms of the standard q, p operators [q, p] = i (see the 
introduction of paper 

Q = q, P = r'^Hsiii{p^^ . (3) 
In what fohows we will restrict our consideration to the first order in /3, that is, we will use the following truncation 

Q = q, P = P+^P^ ■ (4) 

The multidimensional generalization of Eq.([l]) that preserves translation and rotation invariance and introduces a 
finite minimum position uncertainty in all three position variables has the form 0, 0, [l5l | 

[Q\ P'] = I (^--^=J^= S'^ + 2/3P'P^' j , [Q\ = [P\ Pq = . 

The deformed Q, P operators can be represented in terms of the standard q, p operators in the following way 

Its Hilbert space representation (in standard-momentum p representation) has the form 



with the scalar product 



p2</3-i 



To the linear order in /3 Eq.® reads 



= , P,^p,\l + p{pf^ . (6) 
In the classical limit the Eq.® yields the modified dispersion relation 

^p^ +m^ + 2/3p'^ . (7) 

admiting simple physical interpretation that due to quantum-gravitational fluctuations of the background metric, the 
energy e = \/p'^ + rn? acquires the increment /3p^/ -y^p^ -I- m? . 

The corrections to the field theory due to minimum-length deformed quantum mechanics comes in two parts. 

1. Modification at the first quantization level. At the first quantization level the modification amounts to 
the replacement of momentum operator in field Lagrangian with the deformed operator, Eq.®. That is, the modified 
field theory takes the form [13, [3l 

Am = ~ J d'^x^ [^d^<^> + ^P^<P + m^<P^] , (8) 

where P denotes deformed momentum, Eq.® (just to restrict ourselves for simplicity to the leading order in The 
action ([5]) results in the equation of motion 
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9^$ - A<Z> + 2/3 A A<Z> + ni^<P = , 
the plan wave solutions for which can be obtained by the replacement 



exp 1 1 



px — + t 



exp 



px 



/3p4 



(9) 



The equation © explicitly reflects the modified dispersion relation ([7]). In addition to this modification, the theory 
given by the action functional (|8]) involves the cutoff p^ < (that is evident from Eqs.([3l [5])). Once the Lagrangian 
is specified, one then proceeds in the standard way for the second quantization. 

2. Modification at the second quantization level. The second related departure from the usual quantum 
field theory takes place when we are quantizing the field with respect to the deformed quantization [lol , [ll| . Let us 
ignore the modification of the field theory at the first quantization level, Eq.®, and consider second quantization 
with respect to the minimum-length deformed prescription. Assume I is an intrinsic length scale to the problem under 
consideration. After introducing of ladder operators the Hamiltonian of a free field takes the form 



^ = IJ d'fcc^k[fo+(k)6(k)+&(k)fe+(k)] , 
where Wk = (massless field is considered). Defining real variables 



(10) 



the Hamiltonian (jlOp splits into a sum of independent one-dimensional oscillators [ij 



n 



F3 



2Z-1 



■'k^k 



(11) 



each heaving the mass While the appearance of the length scale I in Eq.([TT|) is of no importance in the framework 
of standard quantization (because the energy spectrum of harmonic oscillator does not depend on its mass) , it appears 
explicitly in the energy spectrum when the minimum length deformed quantization is applied [l^ . Namely, assuming 
-Pk, Qk are minimum- length deformed operators and using Eq.Q, the Hamiltonian pTjl to the first order in (3 takes 
the form 



\ J d'k (^ck [fe+(k)6(k) + fe(k)fo+(k)] + [b+ik) - 6(k)]") . (12) 

Thus, at the second quantization level the modification amounts to the replacement of Hamiltonian pUj) with Ea. ([T^ . 
The implementation of modification that arises at the first quantization level into Ea. (ll2l) is now straightforward. 
Namely, one should replace Wk with modified relation that comes out of the modified dispersion relation, Eq. ([7]) . 

The correction term in the Hamiltonian (I12p is controlled by the dimensionless parameter f3/P and therefore the 
validity of our approximation demands the smallness of this parameter. For the correction term involves the length 
scale I, the identification of this scale is a subtle question. In view of Eqs.(Il] [2|) the deviation from the standard 
quantization becomes appreciable at high energies. This fact naturally suggests to identify l^^ with the characteristic 
energy scale of the problem under consideration (for more details see [ll|). For the problem we want to study I is 
naturally set by the Hubble radius H^^. 

Now we are in a position to study the minimum length deformed quantization of a free field on the de Sitter 
background. 



4 



II. MINIMUM-LENGTH DEFORMED QUANTIZATION OF A FREE FIELD IN DE SITTER SPACE 

Let us consider a masslees neutral scalar field in the FLRW background (we will closely follow the paper [isj ) 

To get rid of the overall time-dependent factor, a new field ~ a(r)$ is usually introduced 



A['^] = ^ I drd^x 



<Pdra 



For the Hamiltonian one finds, -/I(x) — dr'P — 'P{dra)/a, 



^ I d^x 



drCl 



Expanding U and 'P in Fourier modes 



77^ + d^Pd^$ + — {n<P + <7>iT) 
a 



where k denotes the comoving momentum, the Hamihonian ()13p takes the form 



•H = I I d^k 



7r(k)7r+(k) + kV(k)(^+(k) + ^ [7r(k)¥.+ (k) + (^(k)^+(k)] 



The quantization conditions 

[<7>(x), i7(y)] = i5{^ - y) , [-Z-Cx), <?(y)] = , [7T(x), i7(y)] = , 
for Fourier amplitudes imply 



(13) 



[^(k), ^(k')] - * <5(k + k') , [^(k), ^(k')] - , [7r(k), 7r(k')] - 



Defining 



6(k) = [u;k^(k) + z^(k)] , 6+(k) = [u;k¥'(-k) - *7r(-k)] , 

where Wk = 'v/k^, one finds 

[6(k), 6+(k')] - 5(k - k') , [6(k), 6(k')] = , [6+(k), 6+(k')] - . 
In terms of &(k), &~''(k) operators the field and momentum operators take the form 



V(2^2L^k 
and the Hamiltonian reduces to 
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c^k [&+(k)6(k) + 6(k)6+(k)] + [5+(-k)6+(k) - 6(k)6(-k)] 



Restricting our attention to the de Sitter space 



J 2 J, 2 2Hti 2 

as — at — e ax , r 



a(r) = - 



the Hamihonian takes the form 



n = - I (fk 



[6+(k)6(k) + fe(k)fe+(k)] + - [6(k)6(-k) - fo+(-k)fe+(k)] 



(14) 



Heisenberg equation of motion 



= i \n, 6(k) = -iwkfc k ^ 

ar r 



has a well known solution 



b{T, k) = 

which for the field operator yields 



6(— oo, k) 



-oo, 



-k), 



(15) 



<?(r, x) = 



v/(27r)32c^k 



6(— oo, k)e 



1- — ) +5+(-oo, k)e*("k^-'^'') + — ) 



(16) 



From the Hamiltonian Eg. (1141) one sees that as r — >■ — cx) the deviation from Minkowski space becomes less and 
less appreciable. This feature allows us to define an unique vacuum state known as the Bunch-Davies vacuum 
6(-cx), k)|0) = 0. 

Introducing real variables (see subsection "Modification at the second quantization level" in Introduction) 



Qk= W^[&(k) + &+(k)] , Pk^ zy^[6+(k)-6(k)] 



the Hamiltonian takes the form 



n = 



(Pk 



k 

2H 



Assuming now Pk, Qk are deformed with respect to the Eq.Q, the Hamiltonian p7)) takes the form 



n 



d^k 



^ + ^j-^(^kP-k+Pkg-k) 



4 P ( 3,3 \ 

^Pk- ^('ZkP^k+Pk'J-kj 



In terms of ladder operators this Hamiltonian reads 



(17) 



(18) 



^ = ^ I d^k 



Wk [^'+(k)6(k) + 6(k)fe+(k)] + - [6(k)&(-k) - 6+(-k)&+(k)] 



12 r 



[6+(k)-&(k)]' 



[6(k) + 6+(k)] [6+(-k) - 6(-k)]^ + [6+(k) - 6(k)]^ [6(-k) + b+{-k) 



(19) 



6 




7 



where now in G(^, k) the foUowing replacement is made 



fit k) ^ e-'--^ 1 



2-key 



The field operator takes the form + /^^i, with given by Eq. p^ and 



<?i(t, k) = 



2wke 



giWk(-r-C) 



1 



2tJkT / 2a;kT \ 2u;kC 



2a;k''' 



1 1 — ^ — — 1 + ^ — - — 1 



G+(e, -k) 

G+(e, k) + 
G(C, -k) 



-ikx 



2a;kCy t^kT 



2tJke y 2a;2rC 2a;kC 



I \ Q-i'^k{-r-(,) jg-iuJk(-r-4) jgiWk(T-^) 



2wk'f y 2u;^Te 2a;ke '^k''' V 2ajkC 



^(2^)32L^k _ 
G(e, k)e^''" -f 
G+(C, k)e-*'' 



X 



(25) 



For the two point spatial correlation function {0\<P{t, 'x)<1>{t, y)|0) (to the first order in /3) one finds 



(0|<?(r, x)<P(t, y)|0) - (0|<Po(t, x)^>o(t, y)|0) + /3(0|<?o(t, x)<?i(t, y)|0) + /3(0|<?i(t, x)<Po(t, y)|0) . (26) 

For both <Po and <Pi are real fields (0|^>o(t-, x)<?i(t, y)|0)* = (0|^>j^(t, x)<?+(r, y)|0) = (0|^>i(t, x)<Po(t-, y)|0). There- 
fore one gets 



(0|<?(r, x)<?>(t, y)|0) = {Q\Mr, x)<Z'o(t, y)|0) + 2/3 5R(0|<Po(t, x)^'i(t, y)|0) 
To the end of inflation, that is r — > 0^, one finds (see Appendix) 



(27) 



\Mr, x)<Pi(r, y)|0)|,^o_ ^ / _ 



sin(kx) sin(ky) cos(kx) cos(ky) 



(28) 



This sort of 5(0) appears also when transition probability is calculated from scattering amplitude in QFT [13] ■ In the 
similar manner we replace it by the three- volume V 



1 



(27r) 



m 



{2n) 



(27r)3 ' 



and set this volume (in view of the problem under consideration) to iJ ^ 

Now returning to the initial field $ = a~^(r)^ and denoting by Tend the end of inflation, one finds 



2/3 5R(0|$o(r, x)$i(t, y 



sin(kx) sin(ky) cos(kx) cos(ky) 



(29) 



The integral in Eq. (j29p diverges in the infrared. To find an appropriate IR cutoff let us estimate the size of pertur- 
bations to the end of inflation 



(0|$2(r, x)|0)|^ 



f dk I3H^ f d^k 



(2^)27 k 256^6^2^^ 



1 sin2(kx) 



12 



(30) 
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that after omitting the sin^(kx) term in Eq. lpO)) takes the form 



m\r, x)|0)|^_ 



dk 



/3H 



.0 



dk 



(31) 



The first integral in Eq. ([3T|) diverges logarithmicahy. The second integral in Eq. ([31]) diverges linearly in the infrared. 
A simple natural way for estimating the scale of an infrared cutoff for the second integral is a follows. The amplitude 
of perturbations in de Sitter space is proportional to the Gibbons-Hawking temperature Tqh = H/2t: J9.]. The 
correction to the Gibbons-Hawking temperature due to minimum-length deformed quantization takes the form 



Tgh — 5> Tgh + I^Tqh 
Therefore, for the corrected perturbation amplitude one obtains 



(32) 



|$2(r, x)|0)| 



[l + (factor of order unity) x 



(33) 



By comparing Ea. ([5T|) with Eg. ([551) one infers that the infrared cutoff should be around ^ 'Tend/ ('^'" somewhat 
smaller - l^^'^r'^d/H as 64 tt^ « 2000). 

It is not out of place to notice that in the case of massive field obeying ^ one finds [l^ 



1$, 



Therefore, in view of Eq. (|32]) one would expect the following correction to Eq. ((33 



(34) 



t; 



GH 



[l -I- (factor of order unity) x /JT^^] 



III. NON-GAUSSIANITIES 



It is clear that the non-linear dependence of the field correction on creation and annihilation operators, Eas. ([^H5|) . 
results in the non-Gaussianity for inflaton perturbations. The non-Gaussianity means that higher order correlation 
functions can not be expressed completely in terms of the two point function [20l |. The lowest order non-Gaussian 
measures arc bispectrum (three point correlation function) and trispectrum (four point correlation function). In the 
case of non-Gaussianity the four point function 



(0|^(t, xi)<P(r, X2)<?(t, X3)<?(r, X4)|0) = (0|<Z'(t, xi)<?(r, X2)|0) (0|<Z'(t, X3)<?(t, X4)|0) + permutations + 

(0|<?(r, xi)<?(r, X2)<P(t, X3)<P(t, X4)|0)cojjj^gctgji , 

contains the connected diagrams (the last term) that vanish in the case of Gaussian perturbations. It is easy to see 
that 

(0|^>(t, Xi)<?(t, X2)<?(t, X3)<?(t, X4) 1 0) connected = ^(0|<?i(t, Xi)<?o(t, X2)<l'o(r, X3)^>o(t, X4)|0) connected + 

permutations with respect to (xi, X2, X3, X4) . 

In our case the field correction (Pi has the following structure, Eqs. 



^>i(x) = B3(x) + C2(x)I?(x) . 

So, using the Wick's theorem in the four point function to the first order in /? we find the following connected 
diagrams 
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(0|B(t, xi)<Z>o(r, X2)|0)(0|B(t, x^-Z-oIt, X3)|0) (0|B(t, xi)<?o(t, X4)|0) , 
(0|C(t, xi)<Po(t, X2)|0)(0|C(t, Xi)<Z'o(t, X3)|0) (0|i^(T, xi)^o(t, X4)|0) . 

(It is worth noticing that the three point function vanishes because in this case the vacuum average is taken over 
an odd number of creation and annihilations operators). Now let us do an explicit calculation. Let us first take the 
term [b'^{£,, — k) — &(^, — k)]^ from Eq. p51) . In Eg. (1^5)) it is accompanied with exp(ikx). We will have the following 
connected diagram 

e^''^^(0| -k) - -k)] <?o(r, X2)|0)(0| [b+{^, -k) - b{^, -k)] <?o(r, X3)|0) x 

(0| [b+it -k) - bit -k)] <^o(t, X4)|0) , (35) 

where the entering vacuum averages can be estimated as 

(oi [6+(e, -k)-6(e, -k)]<?o(T, x2)|o) = 



e 



To estimate the four point function to the end of inflation we may use the following asymptotic expression (see 
Appendix) 

(0| [6+(e, -k) - b{(, -k)] Mr, X2)|0)|^^„_ -Jf^ fT^o ■ (36) 
Thus, using the asymptotic form of Wk(C): Eg. (|46|) . and Eq.(|36|) one finds 



rr4. . ^ikxi 

Wk(C) ^^^^ (0| [b+it -k) - 6(e, -k)] <Po{r, X2)|0) 



(0| [fe+(C, -k) - bit -k)] ^>o(r, X3)|0)(0| [b+it -k) - 6(e, -k)] <Po{r, X4)|0)] 

i/'* /" j3i, sin(kxi) sin(kx2) sin(kx3) sin(kx4) 
96(27r)6r6 J ' 

The other connected diagram contributing to the four point function is (see Eg. pB)) ) 



(37) 



One finds 



e^^^MOl [b+it k) - bit k)] Mr, X2)|0)(0| [b+it k) - bit k)] Mr, X3)|0) x 

(0| [6+(e, -k) + bit -k)] <Po(r, X4)|0) . (38) 



0| [b+it h)-bit k)]<Z'o(r, X2)|0) = 

-il^kC gi(l^kT + kX2) / i \ ( i \ g''("kT-kX2) / j \ 

_ 1 ^ I - g-»"kC ]^ j _ 1 -I I 

2a;kC ^/i2nf2^ \ uj^^r J \ 2wkCy ^/i2nf2U^ \ ^\.r J 

0| [b+it -k) + fe(e, -k)]<^o(r, X4)|0) = 

ii^k? pi(t^kT-kx4) / • \ / -)■ \ „i(i^kT+kx4 



-i ■ 



2a;kC ^(27r)32wk V ^kr/ V 2wk^y ^(27r)32wk V ^kT/ 

So, to the end of inflation one may use the following expressions (see Appendix) 
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(01 k) - k)] Mr. x.)|0)|„„_ ^ J|g= , (39) 

(01 -k) + Mf. -k)] *„(., x,,|0)U._ ^ -±- . (40) 

Again, using the asymptotic form of Wk(0: Eg. (|46l) . and Eqs. (|39l l40|) one finds 



3fJ 



V(27r)32wk 



rr4, . ^ikxi 

1 (0| k) - k)] <|.o(r, X2)|0) 



k) - 6(e, k)] .PoCr, X3)|0)(0| [6+(e, -k) + 6(e, -k)] <Z'o(t, X4)|0)] 



64(27r)6T6 



cos(kxi) sin(kx2) sin(kx3) cos(kx4) 



(41) 



The reason we have taken only real parts in Eqs. ([571HlT) is that after using the complete expression of ^i, see Eq.(|33]), 
only the real parts will survive (and in addition will be multiplied by 2). 

Now returning to the initial field $ = a~^{T)<P and using the infrared cutoff introduced in the previous section, one 
finds the following magnitude for non-Gaussianity to the first order in /3 



\^{Tend, xi)$(re„<j, X2)<i>(Te„d, X3)<i>(re„d, ^4) 1 0) connected 



48(27r) 



6^2 



sin(kxi) sin(kx2) sin(kx3) sin(kx4) + - cos(kxi) sin(kx2) sin(kx3) cos(kx4) 

3 



Estimating order of magnitude for non-Gaussianity (we use r^^^ = iJ^e^^, where N stands for the number of e- 
foldings) 



Non-Gaussianity 



^H^ f d^k _ n(3H^^T^,^^ _ TT^H^e-^^ 



48(27r)6T2„^ J F 48(27r)6 48(27r) 



one infers that, taking into account that for a successful inflation N ~ 70, the non-Gaussianity becomes suppressed 
by the factor e'^^o ^ io-i82^ 



IV. CONCLUDING REMARKS 



In the introduction we tried to clearly outline a scheme for implementation of minimum-length deformed quanti- 
zation in the QFT formalism at the level of second quantization. This scheme is contrasted with the implementation 
of minimum-length deformed prescription in the first quantization picture. Minimum-length deformed quantum me- 
chanics combined with the field theory at the second quantization level reveals qualitatively new features as compared 
to such combination at the first quantization level [lol . [llj , and may teach us something deeper about what is in fact 
implied by this sort of quantization. 

An analytic solution for the field operator is found to the first order in deformation parameter on de Sitter back- 
ground. That is interesting in its own right and can be used for working out of various phenomenological questions. 
The field operator shows up a non-linear dependence on creation and annihilation operators and therefore leads to 
the non-Gaussianity for the perturbation field in the cosmological context. 

Similar solution for the Minkowskian background was found in [ll[. It is instructive to notice that this result 
applied to the black hole radiation yields the correct corrections to the black hole entropy (Tl| . 

Using the solution for field operator, the two point correlation function has been estimated. The correction to the 
two point correlation function shows up an infrared divergence (typical for massless theories). Hence, one is compelled 
to invoke some physical arguments to pin down the infrared cutoff. For this purpose we use the results of paper [ll| 
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that enables one to estimate the order of magnitude for the correction and correspondingly gives an idea about the 
cutoff. The correction to the two point function is of the same order as obtained previously in the framework of 
modified field theory at the first quantization level due to minimum length deformed prescription, see for instance 
most resent paper about it ^9|. 

Non-Gaussianity appears at the lowest level in four point function (the three point function merely vanishes). The 
magnitude of non-Gaussianity is strongly suppressed by the factor e^^^, where N stands for the number of e- foldings. 

It would be interesting to generalize present discussion to the massive field. 

Acknowledgements: Author is indebted to Micheal S. Berger, Jens C. Niemeyer and Zurab K. Silagadze for 
useful comments. 



Appendix 



For notational convenience instead of 6(— oo, k), b^{—oo, k) we simply use fc(k), 6''"(k); besides that we denote 
minimum-length corrected operators by tilde and keep the notations 6(^, k), k) for 



k) = ^ [l ^ j 6(k) + z ^ fo+(-k) , k) = e^--^ + ^ j 5+(k) - . ^-^ b{-k) , (42) 

see Ea. (IT5|) . For calculating the corrections to the correlation function (|26l) . first we have to estimate the following 
vacuum averages (0|fe(k')G(C, k)|0) and {0\b{k')G+ k)|0), where 



G(e, k) = - k) - k)]-^ + ^ (2 -k) - -k)]^ 

+ 3 k) - b{^, k)] ' [6(e, -k) + -k)] + 3 -k) + b+iC, -k)] [b+iC, k) - k)] . (43) 

Splitting Eq.gS]) 

G(^, k) = Gi(e, k) + G2(e, k) , 

Gi(e, k) = - ^ [b+it k) - 6(f, k)]^ + ^ [6+(e, -k) - 6(e, -k)]^ , 

G2(e,k) = :^([6+(e,k)-&(e,k)]'[6(e,-k) + 6+(e,-k)] + [6(^,-k) + 6+(e,-k)] [6+(e,k)-6(e,k)]') , 

one sees that Gi{£_, k) = — Gi(^, k) and G^(f, k) = G2(C, k). So, calculating (0|6(k')G(^, k)|0), we easily estimate 
(0|6(k')G+(^, k)|0) as well 

Non-zero contributions to the (0|6(k') [6+(^, k) - b{^, kW |0) come from the terms containing equal number of 
creation and annihilation operators that cancel each other ^17.] 



(0|&(k')6(-k)6+(-k)fe+(-k)|0) = 2«(5(k' + k)5(0) 



pit^ki pit^kS, 

*7:r^<o|&(k')6+(-k)6(-k)&+(-k)|o) = z5(k' + k)5(o)— — 3 , 

pit^kC / 1 \ pi'^ki I 1 

1 + 2 (0|Kk')W&+(k)&+(-k)|0) = i5{k: + k)5{Q)^—\l 



2wkC V (2wkO / 2wke V (2wkO 

1 + ) (0|6(k')&(-k)6+(-k)6+(k)|0) = - 5{k' - k)(5(0) 



(2wke)' V 2L^ke/ (2wke)' ' 2wke. 

,g-k?h+ ^1 fl + -^Vo|6(k')&+(k)&(k)&+(k)|0) = - 5(k' - k)^(O) e*'^''« ( 1 + ^| fl ' ^ 



(2wkC)V V 2cjke/ V (2wkOV V 2wke/ ' 
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where 5(0) stands for 5(0) = S(\i — k). 

Making replacement k — — k one g 
the above calculations readily determine this term as well. 

Now let us estimate (0|6(k') [b+{^ 
this vacuum average look as follows 



Making replacement k ^ -k one gets (0|6(k') [b+{^, k) - b{^, k)]^ |0) ^ (0|6(k') [b+{^, -k) - b{^, -k)f |0). So, 

Q 

Now let us estimate (0|6(k') k) — &(^, k)]^ — k) + — k)] |0). Terms giving non-zero contributions to 



^ (0|6(k')6+(-k)6(-k)6+(~k)|0) = - (5(k' + k)<5(0) . ' ' 



(0|6(k')&(-k)6+(-k)6+(k)|0) = -iS{k' -k)6{0)- — ^ , 



' ' * (0|6(k')&(k)6+(k)6+(k)|0) = -2z(5(k'-k)(5(0) 



-« ■ 



2wkC V (2wkO / 2wk^ \^ (2a;kO' 

^ h + T^r^ (0|6(k')&+(k)6(k)&+(k)|0) = -*<5(k' - k)<5(0) f- 1 + ——2 



Finally we have to estimate the expression (0|5(k') — k) + — k)] k) — 6(^, k)] |0). One finds that 

it contains the following terms 



(0|6(k')6(-k)6+(-k)6+(-k)|0) = -2S{k' + k)SiO)——^[l-^j 



(2wkO V 2t^ke/ (2wkO 

^ (0|6(k')6+(-k)6(-k)&+(-k)|0) = - d{k' + k)5(0) . ' ^ 



(2a;kO V 2cjkC/ (2a;kO V 2wk^ 



(0|fe(k')&(k)6+(k)6+(-k)|0) -5(k' + k)J(0) 



(2wk^) V 2L^ke/ (2c^kO V 2L^ke 

* (0|6(k')6(-k)6+(-k)6+(k)|0) = -i(5(k'-k)5(0) 



2t^kf y (2cjke) / 2a;ke ^ (2a;kO 



1 + — — (0|&(k')&(k)&+(k)6+(k)|0) = -t2S{k' - k)S{0) 



2wkC V 2a;kC/ 2a;k^ V 2a;kC 

1 + T^-: (0|6(k')&+(k)6(k)6+(k)|0) ^ -i5(k' - k)d{0) ^— 1 



Writing Eg. (1^5)) in a compact form 

T 

<Z>i(r,x) = / f dC[Z^k(e)G(C,k)e^''-+Z^^(C)G+(C,k)e-''-] , (44) 



where 



„ , ., , ,p— '"k(r-C) / • \ iWk(T-4) ,-piWk(T-C) 
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the correlation function (0|^o(''': x)^i(r, y)|0) takes the form 



-iiuji^/T—'k. x) 



(27r)3^u;kii^k' 



[z(0|5(k')Gi(C, k)|0) 5 {Ui^iOe'^^y) + (0|6(k')G2(C, k)|0) (Wk(C)e^''0] 



(45) 



To evaluate this correlation function to the end of inflation, that is, r — > 0^, one can use an asymptotic expression 
for the integral over ^. It is easy to notice that in the limit r — >■ 0~, the main contribution to this integral comes from 
smaller values of that is, ^ — >■ 0. Namely, for a typical term in this integral one finds (n > 1, r < 0) 



< / 



r" " if n is odd , 
if n is even . 



So, to get a leading term when r — )■ , we replace 



T 

(0|6(k')Gi(C, k)|0) 



(0|6(k')G2(e, k)|0) 



4C (2u;kO 



[,5(k' - k) - ,5(k' + k)] (5(0) , 
3 [(5(k' - k) + 5(k' + k)] 6(0) 



After integrating over ^ one gets 



de[*(0|&(k')Gi(^, k)|0)3(Wk(Oe*'' 



(0|6(k')G2(e, k)|0)5R(Wk(Oe*'')] 



48a;2T3 



[S{k' - k) - ,5(k' + k)] sin(ky) + ^^J^^ [S{k' - k) + ^(k' + k)] cos(ky) 



Using the result (|i71) in Eq. (|i5)) one finds 



(46) 



(47) 



||<l'o(r, x)<Pi(r, y 



5(0) 



(27r)3c^k 



48 a;3T4 



sin(kx) sin(ky) + 



2H^ 



64a;3T4 



cos(kx) cos(ky) 



(48) 
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